CONTRACT  TECHNICAL 
REPORT 


"ON  THE  INTERACTION  BETWEEN  DAMAGE  AND 
MOISTURE  SORPTION  IN  CROSS-PLY  LAMINATED 
COMPOSITES" 


By 


Kh.  Abdel-Tawab  and  Y.  Weitsman 


Stotbubow  mxsmmi  n  i 

TBtowd  te  tpa^  leiaoul  | 

Prepared  for:  Office  of  Naval  Research 

Arlington,  Virginia 


Mechanics  Division 
Engineering  Science  Directorate 
Contract  N00014-96-1-0821 


Report  MAES97-1.0-CM 
March  1997 


Mechanical  and  Aerospace  Engineering 
and  Engineering  Science 


THE  UNIVERSITY  OF  TENNESSEE 
Knoxville,  TN  37996-2030 


REPORT  DOCUMENTATION  PAGE 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OeCLASSiFiCATlON/OOWNGRADING  SCHEDULE 


4.  PERFORMING  organization  REPORT  NUMB£R(S) 


64.  NAME  OF  PERFORMING  ORGANIZATION 

DEpt.  of  &  itero.  Ehgr.  &  l 

thiversity  C3f  Itonessee 


6c.  A.OORCSS  (City.  Slau  and  ZIP  Coda} 

307  Fferidus  Ifell 
WxKviUe,  IN  37996-2030 


»4uNAM£  of  FUNOING/SPONSORING 
ORGANIZATION 

Office  of  r^[val  Itesearch 


5c.  ADDRESS  rC/fy.  StaU  ond  ZIP  CcxU) 

Arlingtcxi,  VA  22217 


b.  OFFICE  SYMBOL 
(If  AppUcXibU) 


8to.  OFFICE  SYMBOL 
(If  appUcubU} 


lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION/AVAILABILJTY  OF  REPORT 

LhlimLted 


5.  MONITORING  ORGANIZATION  REPORT  NUM8ER(S) 


74.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  ivfewal  Ifeseardi 


7b.  ADDRESS  {City,  SlaU  and  ZIP  CcxUi 

Office  of  JNfeval  Iteseordi,  Ocxfe  432 
SOON.  &nlryA^;e.,  BOT#!,  Rn.  528' 
Arlingtoi,  m  22217 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


N00014~96-l-0821  (Cm  r 


10.  SOURCE  OF  FUNDING  NOS. 


PjtliiMSgigjgi 


PROGRAM 

PROJECT 

TASK 

YrfORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO. 

12.  personal  AUTHOR(S) 

Khalfid  Abttel-'few^  and  Y.  Jack  Wfeitaran 


134.  TYPE  OF  REPORT 

TfectaLcal 


16.  supplementary  NOTATION 


14.  DATE  OF  REPORT  (Yr.,  Mo..  Dcy) 

15.  PAGE  COUNT 

97-03-03 

9 

COSATI  CODES 


18.  SUBJECT  TERMS  (Continut  on  nti/rrgt  i7  n«c««Mr7  ond  Identify  by  block  number) 


19.  ABSTRACT  (Continut  on  rtotnt  if  nteteeary  and  Identify  by  block  number) 

The  inba3cticn  betweai  itoisture  serptien  and  efenage  in  i±e  fonn  of  transvecse  ccacka  in  cross-^y  laminated 
oenposites  is  studied  assuming  a  linear  Fikian  diffusicn  pcocess.  Finite  elaiHit  cntputaticnal  results/  idnich 
aaxunt  for  the  two-diitaisicnality  cf  tte  pccblan  and  anisotrepy  cf  the  lanunae,  are  cnxpared  against  analytical 
peediGbiens  cf  cne-dimaisicnal  diffusicn.  It  is  ^xwn  that,  for  the  l^-^jp  uniter  oensidaatien,  tie  effect 
cf  transverse  cxad<s  cn  total  moisture  ■weicfit-gain  becares  ncticBable  cnly  iden  the  crack  teeity  ^pxaches  the 
diaracteristic  darter  state  (035). 


20.  OISTRieUTION/AVAILABILlTY  OF  ABSTRACT 

UNCLASSiPtCD/UNuiMlTEO  □  SAME  AS  RPT.  □  OTiC  USERS  □ 

21.  ABSTRACT  SECURITY  CLASSIFICATION 

Dhdassified 

321.  NAME  OF  responsible  INDIVIDUAL 

Dr.  Ifejapakae 

22b.  TELEPHONE  NUMBER 
(Jnetiidt  Area  Cod#; 

(703)  696-4405 

22c.  OFFICE  SYMBOL 

DD  FORM  1473,  83  APR 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


Unclassified 

SECURITY  CLASSIfTcaTION  OF  THIS  PACE 


ON  THE  INTERACTION  BETWEEN  DAMAGE  AND 
MOISTURE  SORPTION  IN  CROSS-PLY  LAMINATED 

COMPOSITES 


Khaled  Abdel- Tawab  and  Y.  Jack  Weitsman 
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ABSTRACT:  The  interaction  between  moisture  sorption  and  damage  in  the  form  of  transverse  cracks 
in  cross-ply  laminated  composites  is  studied  assuming  a  linear  Fikian  diffusion  process.  Finite  element 
computational  results,  which  account  for  the  two-dimensionality  of  the  problem  and  anisotropy  of  the 
laminae,  are  compared  against  analytical  predictions  of  one-dimensional  diffusion.  It  is  shown  that,  for 
the  lay-up  under  consideration,  the  effect  of  transverse  cracks  on  total  moisture  weight-gain  becomes 
noticeable  only  when  the  crack  density  approaches  the  characteristic  damage  state  (CDS). 


INTRODUCTION 

Laminated  composites  are  increasingly  used  in  engineering  applications  that  involve 
exposures  to  ambient  humidities.  Therefore,  for  a  reliable  and  economic  design  of 
structures  made  of  laminated  composites,  it  is  necessary  to  study  moisture  sorption  in 
these  materials  and  the  consequent  effects  on  mechanical  behavior. 

Under  service  conditions,  laminated  composites  may  also  undergo  damage  in  the 
form  of  microcracks,  which  can  affect  the  sorption  process  and  hence  the  moisture 
weight-gain  measurements.  The  objective  of  this  article  is  to  provide  some  guidelines 
for  both  experimental  and  analytical  investigations  of  moisture  content  in  transversely 
cracked  cross-ply  laminated  composites. 

The  effect  of  microcrack  damage  on  the  sorption  process  was  investigated  recently 
[1]  for  tubular  angle-ply  laminated  composites  intended  for  use  in  the  piping  industry. 
In  that  study,  damage  was  associated  with  the  reduction  in  axial  stiffness,  but  was 
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not  physically  observed.  The  main  experimental  finding  was  that  microcrack  damage 
prompted  a  faster  water  absorption  by  the  material  but  retained  the  essential  form 
of  the  absorption  (diffusion)  law  of  the  uncracked  material.  To  account  for  this,  the 
influence  of  damage  was  incorporated  within  an  “effective  time”  that  governed  the 
diffusion  process.  Using  the  concept  of  effective  time,  both  Fickian  and  a  two-phase 
diffusion  [2]  models  were  used  to  fit  the  experimental  data  for  the  moisture  content. 

In  the  present  article,  the  effect  of  transverse  cracks  in  the  90°  ply  group  on  mois¬ 
ture  content  in  cross-ply  laminated  composites  is  investigated  by  considering  moisture 
sorption  to  occur  across  the  transverse  crack  surfaces  as  well  as  the  top  and  bottom 
laminate  surfaces.  Hence  the  sorption  process  is  two-dimensional.  The  present  analysis 
also  accounts  for  the  anisotropy  of  laminated  composites.  For  simplicity,  we  consider 
a  diffusion  process  governed  by  Fick’s  law. 

FORMULATION 

Consider  a  transversely  cracked  [0ni/90n]s  cross- ply  laminate  as  shown  in  Figure  1.  It  is 
assumed  that  the  transverse  cracks  are  uniformly  spaced  at  distance  L  from  each  other. 
The  damage  variable  is  chosen  as  the  transverse  crack  density  a;  which  is  defined  [3]  as 
the  ratio  between  the  thickness  of  the  90°  ply  group  and  the  crack  spacing;  namely 

2  /igo 

~  ~L~ 

Each  ply  group  is  assumed  to  be  transversely  isotropic.  Thus,  for  two-dimensional 
diffusion  in  the  a:ia;2-plane,  Fick’s  law  in  each  0°  ply  group  takes  the  form  [4] 

and  in  the  90°  ply  group  it  takes  the  form 

=  =  (3) 

where  J\  and  J2  are  components  of  the  moisture  flux  vector  in  the  and  X2  direc¬ 
tions,  respectively;  c  is  the  moisture  concentration;  and  and  Dj  are  the  diffusion 
coefficients  in  the  longitudinal  (parallel  to  the  fibers)  and  transverse  (normal  to  the 
fibers)  directions,  respectively. 
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The  diffusion  coefficients  and  Dq-  can  be  determined  from  the  diffusion  coeffi¬ 
cient  of  the  matrix  Dm  and  that  of  the  fibers  Z)y  using  homogeneization  techniques. 
Since  Dy  is  typically  very  small  compared  to  Dm  the  following  approximation  can  be 
used  for  a  fiber  volume  fraction  Vy  <  0.785  [5] 


and 


DL={l-Vq)D, 


DT  =  [\-2.JviJ^)Dm. 


(4) 


(5) 


We  introduce  the  degree  of  anisotropy  7  as  the  ratio  between  Dq,  and  D^ 

It  is  clear  that  7  increases  with  Vy,  and  the  lower  limit  7  =  1  (at  Vy  =  0)  corresponds 
to  an  isotropic  response.' 

Due  to  the  complex  geometry  and  material  anisotropy,  the  finite  element  method  is 
used  to  perform  the  analysis.  In  view  of  symmetry  of  the  problem  it  suffices  to  consider 
only  one  quarter  of  the  domain  between  any  two  adjacent  cracks  as  shown  in  Figure  2. 
It  is  assumed  that  the  transverse  cracks  are  subjected  to  the  same  constant  moisture 
concentration  boundary  condition  (c  =  Co)  as  the  laminate  surfaces.  The  boundary 
conditions  along  with  the  proper  symmerty  conditions  are  shown  in  Figure  2. 

It  should  be  noted  that  the  moisture  concentration  c  remains  finite  at  the  transverse 
crack  tip  at  all  times.  This  follows  from  the  analogy  between  the  present  problem  and 
that  of  heat  conduction  governed  by  Fourier’s  law  [6]  where  it  was  found  that  the 
temperature  remains  finite  at  the  crack  tip. 

The  analysis  was  performed  using  the  ABAQUS  finite  element  program  supple¬ 
mented  with  a  post-processing  subroutine  that  integrates  the  moisture  concentration 
c  to  obtain  the  total  moisture  content  M 

rh  j'L/2 

M{t)  =  c{xi,X2,t)  dxidx2^  (7) 

where  h  =  hQ  +  hgo  is  half  the  laminate  thickness.  In  the  analysis  both  the  transverse 
crack  density  and  the  degree  of  anisotropy  were  varied  by  considering  different  values 
of  u  and  7,  respectively. 


3 


RESULTS  AND  DISCUSSION 


From  the  finite  element  analysis  it  was  found  that  varying  the  degree  of  anisotropy  7 
in  the  range  corresponding  to  0  <  Vy  <  0.785  has  an  insignificant  effect  on  the  total 
moisture  content.  For  all  values  of  the  transverse  crack  density  a;  considered  (Figure  3) 
a  maximum  difference  of  less  than  3%  in  the  moisture  content  is  encountered  over  the 
entire  range  of  7.  Thus,  in  the  remaining  part  of  the  analysis  the  anisotropy  is  neglected 
and  the  response  is  considered  to  be  isotropic  with  a  common  diffusion  coefficient  Dt- 
Different  values  for  the  transverse  crack  density  oo  were  utilized  in  the  finite  element 
runs  and  results  corresponding  to  three  of  these  values  are  shown  in  Figure  3.  The 
results  in  that  figure  are  presented  in  terms  of  a  non-dimensional  time  t*  given  by 


(8) 


where  t  is  the  actual  time.  Also,  the  transient  moisture  content  M  is  normalized  by 
the  steady-state  (maximum)  moisture  content  M{oo)  =  CohL/2.  For  comparison 
purposes,  the  analytical  solution  for  one-dimensional  diffusion  [4]  across  the  laminate 
surfaces  only  (i.e.  in  the  absence  of  transverse  cracks)  is  also  shown  in  Figure  3. 

From  Figure  3a  it  is  clear  that  for  small  values  of  uj  the  effect  of  the  transverse 
cracks  on  moisture  content  is  negligible.  For  intermediate  values  of  to  (Figure  3b)  the 
effect  is  more  pronounced,  but  still  falls  within  the  typical  scatter  band  of  weight-gain 
data  collected  under  one-dimensional  diffusion.  Figure  3c  corresponds  to  a  limiting 
crack  spacing  of  T  =  2/190,  which  approximates  the  so-called  “  characteristic  damage 
state”  (CDS)  of  the  laminate  [7].  In  this  case,  the  effect  of  the  transverse  cracks  is 
quite  considerable  and  can  be  readily  detected  in  an  experimental  program. 

Following  common  practice  in  continuum  damage  mechanics,  damage  can  also  be 
measured  through  the  reduction  in  stiffness.  The  damage  variable  can  then  be  ex¬ 
pressed  as  [8] 


^  Et  ’ 


(9) 


where  Ej  is  the  overall  damaged  axial  modulus  and  Ei  is  the  overall  virgin  axial 
modulus;  both  moduli  measured  in  the  ori-direction.  Suri  and  Perreux  [1]  adopted 
expression  (9)  as  a  measure  for  damage,  and  investigated  moisture  sorption  for  different 
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values  of  a;'  up  to  0.33.  For  transversely  cracked  cross-ply  laminates,  the  ratio  Ei  / 
can  be  determined  using  a  shear  lag  analysis  [3].  Using  material  properties  for  an 
AS4/3501-6  [O/OOaJs  graphite/ epoxy  laminate  [3],  it  is  found  that  the  cases  considered 
in  Figure  3a,  b,  and  c  correspond  to  uj'  «  0.06,0.08,  and  0.17,  respectively.  Larger 
values  for  w'  can  be  realized  by  considering  other  cross-ply  laminates  having  larger 
/jgo  /  Hq  ratios. 

It  is  clear  that  large  variations  in  the  transverse  crack  spacing  correspond  only 
to  small  variations  in  stiffness  reduction.  The  latter  become  even  smaller  for  smaller 
values  of  the  ratio  h^Q  /  Kq.  For  cross-ply  laminates,  it  is  thus  clear  that  expression  (1) 
is  more  appropriate  and  informative  than  (9),  at  least  in  accounting  for  variations  in 
the  transverse  crack  spacing.  On  the  other  hand,  for  angle-ply  laminates  and  other 
complex  laminate  lay-ups  undergoing  distributed  damage,  direct  observations  of  the 
cracks  are  more  difficult  to  perform  and  expression  (9)  could  be  easier  to  adopt  as  a 
damage  measure. 

It  should  be  noted  that  the  idea  of  incorporating  damage  into  “effective  time”  [1], 
does  not  seem  to  be  suitable  for  diffusion  in  cross-ply  laminates  due  to  two  reasons. 
Firstly,  the  presence  of  the  transverse  cracks  renders  the  sorption  process  to  be  two- 
dimensional.  Generally,  a  transient  two-dimensional  boundary  value  problem  cannot 
be  reduced  to  a  transient  one-dimensional  problem  with  an  effective  time.  Secondly,  for 
the  concept  of  “effective  time”  to  be  meaningful,  a  horizontal  shift  function  of  damage 
should  exist  such  that  the  moisture  content  curves  corresponding  to  different  levels  of 
damage  can  be  collapsed  onto  one  master  curve.  This  is  similar  to  the  well-known 
time-temperature  shift  function  used  to  coalesce  creep  curves  at  different  temperatures 
onto  one  master  creep  curve  [e.g.  9].  From  Figure  3,  it  can  be  observed  that  such  a 
shift  function  does  not  exist  except  possibly  during  the  early  stages  of  diffusion  when 
M  is  linearly  proportional  to  \/F.  However,  the  concept  of  effective  time  may  be  useful 
for  diffusion  in  more  complex  laminate  lay-ups. 

Finally,  it  should  be  noted  that  sorption  tests  are  typically  performed  in  two  circum¬ 
stances;  either  under  exposure  to  humid  air  or  with  the  specimens  completely  immersed 
in  fluids.  In  the  latter  circumstance,  if  the  cracks  are  not  drained,  the  weight  of  the 
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liquid  entrapped  in  the  cracks  would  contribute  to  the  measured  weight-gain.  The 
present  analysis  does  not  account  for  such  weight-gain  and  hence  can  be  considered  to 
replicate  moisture  content  data  for  tests  performed  in  humid  air.  In  immersed  envi¬ 
ronments  the  results  presented  herein  remain  valid  provided  that  complete  drainage  of 
the  cracks  is  allowed  prior  to  weight-gain  measurements,  a  process  that  may  take  only 
a  few  minutes  to  complete  in  typical  test  coupons. 
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Figure  3c 


Figure  3.  Normalized  moisture  content  as  a  function  of  non-dimensional  time  for 
different  values  of  transverse  crack  density  lo  {uj  =  2hQofL). 
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